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ABSTRACT 

We report sensitive Spitzer IRS spectroscopy in the 10 — 20 /im region of TW 
Hya, a nearby T Tauri star. The unusual spectral energy distribution of the 
source, that of a "transition object", indicates that the circumstellar disk in the 
system has experienced significant evolution, possibly as a result of planet forma- 
tion. The spectrum we measure is strikingly different from that of other classical 
T Tauri stars reported in the literature, displaying no strong emission features 
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of H 2 0, C2H2, or HCN. The difference indicates that the inner planet formation 
region (< 5AU) of the gaseous disk has evolved physically and/or chemically 
away from the classical T Tauri norm. Nevertheless, TW Hya does show a rich 
spectrum of emission features of atoms (HI, [Nell], and [Nelll]) and molecules 
(H2, OH, CO2, HCO + , and possibly CH3), some of which are also detected in 
classical T Tauri spectra. The properties of the neon emission are consistent with 
an origin for the emission in a disk irradiated by X-rays (with a possible role for 
additional irradiation by stellar EUV). The OH emission we detect, which also 
likely originates in the disk, is hot, arising from energy levels up to 23,000 K above 
ground, and may be produced by the UV photodissociation of water. The HI 
emission is surprisingly strong, with relative strengths that are consistent with 
case B recombination. While the absence of strong molecular emission in the 
10-20 /im region may indicate that the inner region of the gaseous disk has been 
partly cleared by an orbiting giant planet, chemical and/or excitation effects may 
be responsible instead. We discuss these issues and how our results bear on our 
understanding of the evolutionary state of the TW Hya disk. 

Subject headings: (stars:) circumstellar matter — (stars:) planetary systems: 
protoplanetary disks — stars: pre-main sequence — (stars: individual) TW Hya 

1. Introduction 

Spectroscopy with the Spitzer Space Telescope has opened a new window on the gas 
in the planet formation region of circumstellar disks (< 5AU). It has revealed the presence 
of organic molecules in this inner region of the disk (Lahuis et al. 2006b; Carr & Najita 
2008) and shown that water emission from disks is both more commonly occuring and ex- 
tends to larger radii than previously known from ground-based observations (Carr & Najita 
2008; Salyk et al. 2008). Spitzer spectroscopy has also demonstrated that [Nell], which was 
predicted to arise from irradiated disk atmospheres (Glassgold et al. 2007), is commonly 
present in spectra of T Tauri stars (Pascucci et al. 2007; Lahuis et al. 2007; Espaillat et al. 
2007; Guedel et al. 2009). These diagnostics offer the opportunity to study the evolution of 
gaseous disks and to obtain new insights into the processes by which disks evolve. 

Transition objects are an interesting class of sources to study in this context. Their 
spectral energy distributions (SEDs) are interpreted as arising from an optically thin inner 
disk (within an opacity "hole" of radius i?hoie) that is surrounded by an optically thick outer 
disk (beyond -Rhoic)- Such an SED can arise in several ways, as a consequence of grain 
growth (e.g., Strom et al. 1989; Dullemond & Dominik 2005), giant planet formation (e.g., 
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Skrutskie et al. 1990; Marsh & Mahoney 1992; Lubow et al. 1999; Bryden et al. 1999), or 
disk photoevaporation (Clarke et al. 2001; Alexander et al. 2006). Transition objects might 
arise in any of these ways, and the evolutionary path for any given T Tauri star, which may 
or may not involve a stint as a transition object, may depend on initial conditions such as 
the initial disk mass. 

The above processes do make different predictions for the stellar accretion rates and 
disk masses under which a transitional SED will arise (Najita et al. 2007a; Alexander & 
Armitage 2007). As a result, stellar accretion rates and disk masses have been employed 
previously in diagnosing the nature of transition objects (Najita et al. 2007a; Cieza et al. 
2008; Alexander & Armitage 2009). 

Because these processes also make different predictions for the structure of the gaseous 
disk, studying the gaseous disks of transition objects is potentially another tool that can be 
used to distinguish among these possibilities (Najita et al. 2007a, 2008). This approach, the 
one taken here, may therefore offer insights into how planets form and disks dissipate. 

At a distance of 51 pc (Mamajek 2005), TW Hya is the nearest T Tauri star and a well- 
known transition object. Its SED indicates that the region of the disk within i?hoie ~ 4AU 
of the star is optically thin in the continuum (Calvet et al. 2002; Uchida et al. 2004). Higher 
angular resolution interferometric results at 7mm (Hughes et al. 2007) and at 10 /im (Ratzka 
et al. 2007) further confirm a deficit of continuum emission from the central region of the 
disk. TW Hya is also an unusual T Tauri star in other respects. Despite its advanced 
age (~ 10 Myr), it has a massive disk (Wilner et al. 2000; Calvet et al. 2002) and is still 
accreting (Muzerolle et al. 2000; Johns-Krull et al. 2000; Herczeg et al. 2004; Alencar & 
Batalha 2002). TW Hya is also a source of energetic UV emission (Herczeg et al. 2004; 
Bergin et al. 2004) and emits strongly in X-rays (Kastner et al. 1999; Stelzer & Schmitt 
2004; Robrade & Schmitt 2006). 

Spectroscopy of TW Hya taken with the Spitzer Infrared Spectrograph (IRS) has been 
reported previously, with an emphasis on the shape of the mid-infrared continuum and what 
it says about the radial distribution of the dust opacity in the disk (Uchida et al. 2004). Based 
on a reanalysis of the same data, Ratzka et al. (2007) reported the detection of emission lines 
of HI 6-5, HI 7-6, and [Nell]. The [Nell] emission from TW Hya has been studied from the 
ground at high spectral resolution by Herczeg et al. (2007) and Pascucci & Sterzik (2009). 

Here we report further IRS spectroscopy of TW Hya obtained at higher signal-to-noise 
in two later epochs. The spectra show a rich suite of emission features, both atomic and 
molecular, some of which have been detected before in spectra of classical T Tauri stars. 
We discuss the physical processes that may be responsible for the emission (section 3). The 
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spectra are also strikingly different from the spectra of classical T Tauri stars, in that they 
lack strong emission features of H 2 0, C 2 H 2 , and HCN. This suggests that the inner gaseous 
disk has evolved physically and/or chemically away from the classical T Tauri norm. We 
discuss how these differences bear on our understanding of the evolutionary state of the disk 
(section 4.1). We also discuss where in the system the detected emission features may arise 
(section 4.2). 



TW Hya was observed twice with the Short-High (SH) module of the IRS (Houck et al. 
2004) on board the Spitzer Space Telescope (Werner et al. 2004) as part of the program GO 
30300. A first epoch spectrum (hereinafter SHI) was acquired in 2006. A second, longer 
integration spectrum (hereinafter SH2) was acquired in 2008. An earlier, shorter integration 
spectrum that pre-dates both of the spectra from our GO program was acquired in 2004 
as part of the IRS GTO program (hereinafter GTO). The GTO data have been presented 
previously by Uchida et al. (2004) and Ratzka et al. (2007). Table 1 gives a log of the 
observations. 

All three of the data sets were reduced using the techniques described in Carr & Najita 
(2008). The GTO dataset lacks off-target frames and consists of just one exposure in each of 
the nod positions. In this case, the identification of bad pixels relied solely on the inspection 
of the single images, using a comparison of the two nod positions as an aid. The spectra 
were calibrated by creating a spectral response function from a suite of archival standard 
star spectra reduced in the same way as the target spectra. The spectral response function 
was customized for each nod position of each data set by by combining those standard 
star spectra that minimzed the fringing in order 20. Subsequent application of IRSFRINGE 
(Lahuis et al. 2006a) found little residual fringing. Table 2 provides estimates of the resulting 
noise per pixel and the continuum level at the wavelengths of selected spectral features in 
the 10 — 20 fim range. The noise in each spectral region was estimated from the difference 
of the spectra taken in the "A" and "B" beam positions at each epoch. 



2. 



Observations and Data Reduction 



3. 



Results 



3.1. 



Continuum and Emission Lines 



The shape of the mid- infrared continuum is similar in all three spectra (Figure 1). 
The continuum in the GTO spectrum is well approximated by 1.1 times the flux of the 
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continuum in SH2. The SHI continuum is approximately 1.3 times the SH2 continuum, 
although it is shallower in slope than SH2 beyond ~ 15 fim. These differences may arise in 
part from pointing errors (< 20% photometric uncertainity in SH for the nominal Spitzer 
pointing accuracy; see Spitzer Observer's Manual, version 8.0) or true variability in the T 
Tauri continuum (e.g., Muzerolle et al. 2009) or a combination of these. The [Nell] line at 
12.81 fim and HI lines at 11.31 fim and 12.37/xm are strongly in emission at each epoch. 

Notably, in none of the three epochs does TW Hya show strong mid-infrared molecular 
emission like that seen in Spitzer spectra of classical T Tauri stars that have been reported 
in the literature (Carr & Najita 2008; Salyk et al. 2008; see also Pascucci et al. 2009). The 
[Nell] line luminosity is a convenient reference, since TW Hya and AA Tau (Carr & Najita 
2008) have similar [Nell] line luminosities, within a factor of 2. In AA Tau, the molecular 
emission is stronger than the [Nell] emission, whereas in TW Hya, the [Nell] emission is, by 
far, one of the brightest lines in the spectrum. 

Nevertheless, TW Hya does show a rich spectrum of weaker emission features. To 
illustrate this, we show in Figure 2 the average of SHI and SH2 after subtracting a high order 
polynomial fit to the continuum. The features that are seen correspond well in wavelength 
with known atomic lines (HI, Nell, Nelll) and other molecular features (H 2 , OH, C0 2 , 
HCO + , and possibly CH3). The transitions that dominate the detected OH features are 
listed in Table 3. The CO2 band has been detected both in emission (Carr & Najita 2008; 
Salyk et al. 2008) and in absorption (e.g., Lahuis et al. 2006b) from T Tauri stars. 

We identify the feature at 12.06 fim as the z/ 2 fundamental Q-branch of HCO + (Davies 
& Rothwell 1984). The detection of HCO + emission is interesting, and to our knowledge 
this band has not been reported previously in any young stellar object. There is also a 
weak feature observed in the SHI spectrum at 16.48 fim that may be the v 2 Q-branch of the 
fundamental band of CH 3 (Bezard et al. 1999; Feuchtgruber et al. 2000). In addition, the 
[Nelll] line is marginally detected in the SH2 spectrum. Figure 3 gives an expanded view 
of the region around the weak features [Nelll] and CH 3 in the epoch in which they were 
detected (SH2 and SHI, respectively). These features are more apparent in the individual 
spectra than in the average of SHI and SH2 (Figure 2). 

Since fitting and subtracting a high order polynomial can affect the strength of weak 
features (e.g., if it incorrectly fits the local continuum), we did not measure line strengths 
from the continuum-subtracted spectrum. We instead used the continuum-subtracted spec- 
trum to identify the approximate wavelengths of features to be fit, and we then measured 
feature strengths from the original spectrum. We made a least-squares fit to each feature 
using a model of a Gaussian and a quadratic approximation to the local continuum. A 
Marquardt routine was used to find the best fit (e.g., Bevington 1969). Estimated errors on 
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the line fluxes and line widths were found by determining the range of values enclosed by 
the 68% confidence interval (a "l-cr" estimate) that is appropriate for the 6 parameters in 
the fit (e.g., Lampton et al. 1976). 

Tables 4 and 5 report the fluxes, Gaussian widths, and equivalent widths of the features 
in the SHI and SH2 spectra, respectively. Table 6 reports the fluxes and equivalent widths 
measured from fits to the features in the GTO spectrum. Feature widths are not reported 
for the GTO spectrum because they are not well constrained given the lower signal-to-noise 
of that spectrum. The fluxes of the features measured in the GTO spectrum (HI 9-7, HI 
7-6, [Nell], H 2 , and HCO + ) are consistent within the errors with the fluxes measured in SHI 
and SH2. 



3.2. Properties and Variability of Emission Features 

For the lines that appear to be unblended (i.e., the HI lines 12-8 at 10.5 /im, 9-7 at 
11.31 /im, 7-6 at 12.37 /im, 13-9 at 14.18 /im, 8-7 at 19.07 /im, [Nell] at 12.81 /im, and the 
H2 SI line), the line center positions are typically within ±50kms _1 of their rest values. As 
shown in Figure 4, the widths of these lines are, for the most part, consistent with being 
spectrally unresolved. 

We detect a series of OH features that are each comprised of multiple pure rotational 
transitions of OH. Table 3 gives the approximate wavelengths of these OH groups and the 
transitions that dominate the flux we observe. Most of the emission arises from transitions 
within the v =0 ground vibrational state, from both the 2 n 3 / 2 and 2 Hi/2 rotational ladders, 
with J=14.5 to J=29.5. Figure 3 of Tappe et al. (2008) gives a graphical illustration of the 
range of upper energy levels probed by these transitions. Some low lying v=l transitions 
may also contribute by enhancing the v =0 emission groups when there is a wavelength 
coincidence (e.g., at 14.65 /im and 15.30 /im). The widths of the OH features are comparable 
to or broader than the nominal spectral resolution of IRS; e.g., the features at 14.08 /xm, 
14.65 /xm, 15.30 /im, 16.02 /mi, and 17.78 /im are spectrally resolved (Figure 4). The increase 
in the feature widths with wavelength arises because the spread in the wavelengths of the 
lines that contribute to each feature is progressively larger at longer wavelengths. 

Figure 5 shows the ratio of the line fluxes, equivalent widths, and Gaussian line widths 
in SHI compared to SH2. The flux of the [Nell] line is similar in the two epochs (within 
~ 10%), but its equivalent width is lower (by ~ 30%) in SHI than in SH2 because of 
the brighter continuum in SHI (Figure 1). In comparison, the HI lines at 11.31 /xm and 
12.37 /xm are stronger in both flux (by a factor of ~ 2) and equivalent width (by ~ 30%) in 
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SH1 compared to SH2. The variation of the HI lines suggests true variability. Neither the 
12.37/im HI line nor the [Nell] line are spatially extended with respect to the continuum at 
the spatial resolution of Spitzer, hence, while a small pointing error could alter the measured 
line flux, it is unlikely to alter the line equivalent width. The 12.37 /xm line is brighter than 
the [Nell] line in SHI and fainter than [Nell] in SH2, also indicating true variability. The 
H2 lines have similar fluxes in both epochs, as does HCO + . 

Some features are apparent in one epoch but not the other. We find in SH2, but not 
SHI, a weak feature at the wavelength of the [Nelll] 15.55 /im line. The [Nelll] line may 
be more apparent in SH2 because of the higher signal-to-noise of the SH2 spectrum. We 
find in SHI, but not SH2, weak features at the wavelength of the CH3 band (at 16.48 /im) 
and the HI 14-9 line (at 12.58 /im). Finding the latter feature in SHI but not the higher 
signal-to-noise SH2 is consistent with the brighter emission observed for the dominant HI 
lines (7-6 and 9-7) in SHI than in SH2 (Figure 1). 



3.3. HI Lines 

The flux of the HI 7-6 line, as an average of the flux in SHI (2.3 x 10 28 ergs -1 ) and SH2 
(1.3 x 10 28 ergs _1 ), is in good agreement with the HI 7-6 flux of 1.7 x 10 28 ergs _1 reported 
by Ratzka et al. (2007). 

The relative fluxes of the HI lines in SHI and SH2 can be accounted for roughly with 
case B recombination. Figures 6 and 7 compare the HI fluxes measured in SHI and SH2, 
respectively, with the relative fluxes for case B at a temperature of T = 10 4 K and an electron 
density of n e = 10 4 cm -3 (Storey & Hummer 1995). The case B relative fluxes are similar 
(i.e., within our measurement errors) over a broad range in temperature (1, 000 — 20, 000 K) 
and electron density (n e = 10 3 — 10 8 cm -3 ). As a result, the relative fluxes we measure do 
not constrain the temperature and electron density of the emitting gas. 

In comparison, greater variation is expected in the relative line fluxes of lower lying HI 
lines over the same range in density and temperature. As described by Bary et al. (2008), 
the Paschen (Pa/3, Pa7, Pa5, and 8 through 14) and Brackett (Br7 and 10 through 16) lines 
of T Tauri stars that originate from the same range of upper energy levels also show line 
ratios consistent with case B recombination. They interpret their line ratios as favoring low 
temperature (< 2000 K) and high electron density (10 9 cm -3 < n c < 10 10 cm -3 ). 

The mid-infrared (MIR) lines we observe are predicted to have inverted populations at 
such high densities (Hummer & Storey 1987). Transitions between high lying levels typically 
become inverted at lower densities than those between low lying levels. So, for example, the 
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HI 9-7 and 8-7 transitions have a population inversion at densities above n e ~ 10 8 cm" 3 , 
whereas the HI 7-6 transition is inverted at higher densities (n e > 10 8 — 10 9 cm -3 ) and 
the HI 13-9 transition at lower densities (n e > 10 7 cm" 3 ). The fact that the relative fluxes 
we observe in TW Hya are consistent with case B recombination suggests that we are not 
seeing maser emission in this source, as the predicted amplification factors can differ widely 
between transitions. Similarly, MIR HI transitions were not found to be masing in a study 
of MWC349 (Smith et al. 1997), a source that shows maser emission in longer wavelength 
HI lines. 

We can estimate the number of ionizing photons that are needed to produce the observed 
HI emission under the assumption of ionization equilibrium, i.e., that the number of ioniza- 
tions per second equals the number of recombinations per second. Since the recombination 
line luminosity L u[ = n e n p e u iV (where V is the optically thin emitting volume, n e is the 
electron density, and e u i is the emissivity of the transition in ergs -1 cm 3 ), the total recombi- 
nation rate aBfi e n p V can be expressed as as L u i/e u i. For the HI 9-7 transition at 11.31 fim in 
SH2, L ui = 4.5 x lO^ergs" 1 . From Hummer & Storey (1987), e ul = 3.7 x KT^ergs" 1 cm 3 , 
and as = 2.6 x 10 _13 cm 3 s _1 for T e = 10 4 K and n e = 10 4 cm -3 . For these values, the total 
recombination rate is 3.2 x 10 42 s _1 . 

In comparison, the ionizing photon flux of T Tauri stars is often estimated as 10 41 s" 1 
(e.g., Hollenbach et al. 2000), although Alexander et al. (2005) derived ionizing photon fluxes 
of 10 41 — 10 44 s _1 for various T Tauri stars based on CIV A1549 line luminosities. Looking 
at the situation for TW Hya itself, Herczeg (2007) estimated an ionizing photon flux of 

5 x 10 41 s" 1 from the accretion shock onto the star based on the observed X-ray and FUV 
line fluxes and the assumption of collisional ionization, with a factor of 5 uncertainty in the 
quoted value. The number of ionizing photons that we infer is at the upper end of this range, 
so TW Hya may emit enough ionizing photons to explain the mid-infrared HI line fluxes that 
we observe. However, if the EUV photons are sufficient to produce the mid-infrared HI line 
fluxes, they are likely to produce more [Nell] or [Nelll] emission than is observed (Hollenbach 

6 Gorti 2009; see also below). 

We can also estimate the mass of the emitting (ionized) gas, under the assumption of 
optically thin emission and constant density. The mass of ionized hydrogen is van n p V = 
(ran / n e )(L u i / e u i) = 1.6 x 10 31 gcm _3 /n e . This corresponds to a total number of hydrogen 
atoms of 10 55 cm~ 3 /n c . If we require n e < 10 8 cm" 3 to avoid inverted populations, the 
emitting mass of ionized gas is > 10 23 g. In section 4, we consider the region of the YSO 
system in which such emission might arise. 
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3.4. Neon Lines 

The range of [Nell] equivalent widths we measure is consistent with the equivalent width 
reported for this line from ground-based high resolution spectroscopy. We measure equivalent 
widths of 63±10A, 44±3A, and 65±2A in the GTO, SHI and SH2 spectra, respectively. In 
comparison, Herczeg et al. (2007) measured at high spectral resolution an equivalent width of 
62 ± 11 A for the [Nell] line using MICHELLE on Gemini. TW Hya was also studied at high 
spectral resolution by Pascucci & Sterzik (2009) using VISIR on the VLT. They measured a 
[Nell] equivalent width of 4lA and a flux of 3.8 x 10~ 14 ergs _1 cm~ 2 , within 65-75% of the 
flux of 5.9 x 10~ 14 erg s _1 cm -2 , 5.0 x 10~ 14 erg s _1 cm -2 and 5.6 x 10 -14 erg s _1 cm~ 2 measured 
for the GTO, SHI, and SH2 spectra, respectively. 

Time variability, pointing errors, flux calibration issues, and/or the narrower slit of the 
ground-based observations may account for any differences between the IRS and ground- 
based results. In particular, the 0.4" slit of the VISIR observations translates into a physical 
size scale of ±10 AUat the 50 pc distance of TW Hya. As noted by Pascucci & Sterzik (2009), 
if the [Nell] emission arises from a disk, the narrow slit may exclude a reasonable fraction of 
the emission, since Glassgold et al. (2007) have shown that the [Nell] emitting region of the 
disk may extend out to a radius of ~ 20 AU. We further note here that in the calculations 
of Meijerink et al. (2008) for their fiducial T Tauri disk, ~ 75% of the [Nell] emission arises 
from within 10 AU. If there is no substantial gaseous disk within the optically thin region 
of the TW Hya disk (e.g., if the disk is dynamically cleared by an orbiting planet), then it 
would be more appropriate to look at the emission arising from 4-10 AU. The fraction of the 
[Nell] emission that arises from this region of the disk is ~ 65%, comparable to the fraction 
of the IRS [Nell] flux that is recovered in the VISIR data. 

The MICHELLE observations of Herczeg et al. (2007) were also made with a 0.4" slit, 
but in poor seeing. So slit losses may be an issue there, but the observations are likely to 
be less sensitive to a difference in the spatial extent of the [Nell] emission relative to the 
MIR continuum. That is, the equivalent width of the observations could be similar to that 
measured by Spitzer, as observed. Near-simultaneous narrow- and wide-slit ground-based 
[Nell] observations of TW Hya, taken in good seeing, would be useful in constraining the 
spatial extent of the emission. 

While [Nell] is commonly detected from T Tauri stars (Pascucci et al. 2007; Espaillat 
et al. 2007; Lahuis et al. 2007; Ratzka et al. 2007; Guedel et al. 2008, 2009; Flaccomio et 
al. 2009), [Nelll] is rarely detected, with upper limits typically an order of magnitude lower 
(e.g., Lahuis et al. 2007). In addition to our possible detection of [Nelll] from TW Hya, 
the T Tauri stars Sz 102 and WL5 also show [Nelll] emission. Sz 102 has neon fluxes and 
a [NeIII]/[NeII] ratio (~ 0.06; Lahuis et al. 2007) that are similar to the values reported 
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here for TW Hya ([NeIII]/[NeII]~ 0.045). In contrast to TW Hya, Sz 102 drives a jet (e.g., 
Coffey et al. 2004, where the source is known as TH 28), so its neon emission may have 
contributions from a both a disk and a jet. [Nelll] is also detected from WL5, where the 
[NeIII]/[NeII] flux ratio is much higher, ~ 0.25 (Flaccomio et al. 2009). 

The [Nell] and [Nelll] line fluxes that we detect from TW Hya are in good agreement 
with predictions for gaseous disk atmospheres that are irradiated by stellar X-rays. For 
an X-ray luminosity of 2 x 10 30 ergs _1 , which is roughly the X-ray luminosity of TW Hya 
(Kastner et al. 1999; Robrade & Schmitt 2006), the models of Glassgold et al. predict a [Nell] 
luminosity of 1.4 x 10 28 erg s _1 and a [Nelll] luminosity of 1.6 x 10 27 erg s -1 . In comparison, we 
measure [Nell] and [Nelll] luminosities of 1.6 x 10 28 ergs _1 and 0.7 x 10 27 ergs _1 in the SH2 
spectrum. In the Glassgold models, the [NeIII]/[NeII] ratio is small because the emission 
arises in a region with a significant abundance of neutral hydrogen atoms (the ionization 
fraction x e ~ 0.001). As a result, there is efficient charge exchange between Nelll and HI 
and a consequent reduction in the Nelll abundance. 

EUV photons reaching the disk would enhance the neon emission above that produced 
by X-ray irradiation. For an EUV irradiated region, which is typically fully ionized (x e ~ 1), 
Hollenbach & Gorti (2009) find that a small [NeIII]/[NeII] ratio can be produced with a 
soft EUV spectrum, e.g., a blackbody with a temperature of 30,000 K. A harder spectrum, 
e.g., with constant uL u , would produce [NeIII]/[NeII] ~ 1, in conflict with the observations. 
The EUV spectrum of TW Hya has not been measured directly, although Herczeg (2007) 
has estimated a significant EUV flux from the accretion column based on the observed 
FUV and soft X-ray spectrum. Herczeg did not comment on the spectral shape of the 
predicted EUV emission. The impact of the emitted EUV on the disk may be limited by 
foreground HI absorption (e.g., in a magnetosphere; Alexander et al. 2004). Reviewing the 
literature, Herczeg (2007) estimated that only a small fraction of the EUV (equivalent to 
10 39 photons s" 1 ) reaches the disk. So the EUV contribution to the neon emission from the 
disk may be limited. 



3.5. Molecules 

The OH features we detect originate from very high energy levels (upper energy lev- 
els 6,000 K to 23,000 K above ground) and have large transition probabilities (A-values of 
20 — 70s -1 ). OH emission from other young stellar disks has been reported previously for 
lower rotational lines in both the ground vibrational state (Carr & Najita 2008) and the 
fundamental ro- vibrational band (Mandell et al. 2008; Salyk et al. 2008). The rotational 
temperatures derived for the emission from these lower rotational levels are consistent with 



- 11 - 



thermal OH emission at 500-1000 K. The much higher energy rotational states observed in 
emission in the SH spectrum of TW Hya are more similar to the superthermal OH emission 
detected by Tappe et al. (2008) in the Spitzer IRS spectrum of young stellar outflow HH 
211. 

One possible interpretation of the highly excited OH emission detected from TW Hya 
is that it arises from gas that is hot and dense enough to collisionally populate the upper 
levels. The excitation conditions for such thermal OH emission would exceed even those for 
the CO fundamental emission observed from TW Hya (upper energy levels 3,000 to 6,000 K 
and A-values ~ 10s _1 ), which arises from gas at ~ 1000 K in the inner 1 AU of the disk 
(Salyk et al. 2007). The ro- vibrational transitions of CO are generally considered to be a 
good tracer of hot and dense molecular gas in the circumstellar environment, and the lack 
of a CO component corresponding to the hotter OH emission is curious. 

The lack of H2O emission, along with a highly excited OH emission spectrum, is a 
striking difference between TW Hya and normal classical T Tauri stars. However, close 
examination of the IRS SH spectra of classical T Tauri stars does reveal the presence of these 
same OH emission features (Carr & Najita, in preparation), with rotational temperatures 
~ 4000 K. This emission was not previously recognized because of the complexity of the 
dense H2O emission in which the OH emission is embedded. Hence, the character of the 
OH emission in TW Hya and that of a classical T Tauri star such as AA Tau may be more 
similar than would appear at first glance. 

A likely possibility is that such highly excited OH emission arises from the photodis- 
sociation of H2O, which is also the interpretation given by Tappe et al. (2008) for the OH 
emission from HH 211. The photodissociation of H2O by UV photons in the second ab- 
sorption band of water (A = 1200 — 1400A) will produce OH molecules in highly excited 
rotational levels but predominantly in the ground vibrational state (Harich et al. 2000). We 
will present a detailed analysis of the OH spectrum of TW Hya in a following paper. 

The H2 emission that we observe is much weaker than either the [Nell] or HI 7-6 
emission from TW Hya. The line luminosities are in the range found for other accreting 
young stars. In the ground-based (IRTF/TEXES) H 2 survey of Bitner et al. (2008), the H 2 
S(2) line luminosity of a low mass star such as GV Tau N is ~ 3 x 10~ 6 L©, although upper 
limits on the S(2) line luminosity for some T Tauri stars are much lower (e.g., ~ 0.4 x 10 _6 L© 
for LkCa 15). In comparison, the line luminosities of the S(l) and S(2) emission that are 
detected from TW Hya in SH2 are 0.8 x 10~ 6 L© and 0.5 x 1O~ 6 L , respectively. Since the 
S(2) line of TW Hya is blended with an OH feature, we assumed an OH contribution that 
is similar to that of neighboring OH features in estimating the flux of the S(2) line. The 
observed S(l) luminosity is similar to the value of 0.2 x 10~ 6 L Q that is predicted for TW 
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Hya by Gorti & Hollenbach (2008). 

Interestingly, the H 2 emission from TW Hya is stronger than the upper limits on the 
H 2 flux that Bitner et al. (2008) measured for TW Hya. Whereas we measure line fluxes of 
1.3 x 10- 14 ergs" 1 cm" 2 and 0.6 x 10~ 14 ergs" 1 cm" 2 for the S(l) and S(2) lines (5-a and 3-a 
detections, respectively, in SH2), Bitner et al. (2008) provide 3-a upper limits for the S(l) 
and S(2) lines of 0.7 x 10 _14 ergs _1 cm~ 2 and 0.6 x 10 _14 ergs _1 cm~ 2 , respectively. Perhaps 
the narrower slit (0.81" and 0.54" for the SI and S2 lines, respectively) of the TEXES 
measurements plays a role here as well. The line widths that are measured for H 2 emission 
from young stars, when it is detected by TEXES in more distant sources (~ 140 pc), are 
~ lOkms^ 1 , which is consistent with much of the emission arising from the 10-50 AU region 
of the disk. As discussed in the context of the [Nell] emission from TW Hya (section 3.4), 
much of this region of the disk is excluded by the narrow slit in a nearby source such as 
TW Hya. As in the case of [Nell], near-simultaneous narrow and wide slit ground-based 
observations of the H 2 emission from TW Hya, taken in good seeing, would be useful in 
constraining the spatial extent of the emission. 

We also detect molecules other than OH and H 2 in the TW Hya spectrum (C0 2 , HCO + , 
and possibly CH 3 ). The properties of these emission features will be discussed in greater 
detail in a future publication. 

4. Discussion 

4.1. Evolutionary State of the TW Hya Disk 

As a well-studied, nearby transition object, TW Hya has been a touchstone for under- 
standing the origin and evolutionary state of such systems. Transition objects have unusual 
spectral energy distributions characterized by a deficit of infrared excess at short wavelengths, 
a trait that is interpreted as indicating that the disk is optically thin in the continuum within 
a given radius i4oie (Strom et al. 1989). Analyses of spectral energy distributions of transi- 
tion objects find -Rhoie values of a few AU to ~ 50 AU in some cases (Espaillat et al. 2007). 
In comparison, the spectral energy distribution of TW Hya indicates that its disk is optically 
thin within ~ 4 AU of the star (Calvet et al. 2002; Uchida et al. 2004). An inner hole in the 
dust continuum emission is detected at 7mm and at 10 /im (Hughes et al. 2007; Ratzka et 
al. 2007). 

Such a dust distribution could be interpreted as a result of grain growth (Strom et al. 
1989; Dullemond & Dominik 2005), photoevaporation (Clarke et al. 2001; Alexander et al. 
2006; Gorti & Hollenbach 2009; Owen et al. 2010), or a gap opened by an orbiting giant 
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planet (e.g., Marsh & Mahoney 1992; Calvet et al. 2002; Rice et al. 2003; Quillen et al. 
2004; Calvet et al. 2005). These scenarios could potentially be distinguished by examining 
properties other than the SED, since the scenarios make different predictions for the stellar 
accretion rate, disk mass, and radial distribution of gas in the disk (Najita et al. 2007a; 
Alexander & Armitage 2007). 

The first two diagnostics (stellar accretion rates and disk masses) were used by Najita et 
al. (2007a) to probe the nature of transition objects in the Taurus star forming region. They 
found that Taurus transition objects have higher than average disk masses as well as stellar 
accretion rates that are ~ 10 times lower than non-transition objects. These properties are 
roughly consistent with the predictions of theories of giant planet formation (e.g., Lubow 
et al. 1999; Lubow & D'Angelo 2006). The high disk mass of TW Hya (> O.O6M ; Calvet 
et al. 2002) and its comparatively low stellar accretion rate (~ lO~ 9 M yr _1 ; e.g., Herczeg 
et al. 2004; Muzerolle et al. 2000; Alencar & Basri 2000) place it in a similar region of the 
M-M disk plane as the Taurus transition objects. 

Studies of line emission from the gaseous component of transition disks, like that pre- 
sented here, offer the opportunity to complement studies of stellar accretion rates and disk 
mass, by probing the radial distribution of gas in the disk and therefore the evolutionary 
state of the system. As described by Najita et al. (2007a, 2008): 

(1) In the grain growth and planetesimal formation scenario, the inner disk is rendered 
optically thin in the continuum, but the gaseous component is unaltered and would fill 
the region within -Rhoie- Emission from gas everywhere within -Rhoic would produce bright 
emission because of the lack of continuum emission from same region of the disk. 

(2) If a planet has formed with a mass sufficient to open a gap (~ lMj), gas will be cleared 
in the vicinity of its orbit, but gap-crossing streams, from the outer disk to the planet, and 
from the planet to the inner disk, can allow continued accretion onto both the planet and 
the star, the latter via the replenishment of the inner disk within -Ri nne r < -Rhoic (e.g., Lubow 
et al. 1999; Kley 1999; Bryden et al. 1999; D'Angelo et al. 2003; Bate et al. 2003; Lubow 
& D'Angelo 2006). While the inner disk (r < R mnC r < Rhoic) might then produce gaseous 
emission lines, the low surface filling factor of gas in the region of the gap would produce 
weak to negligible emission because of the small projected emitting area of the accretion 
streams. 

(3) In the case of photoevaporation, gas in the region of the disk photoevaporation radius 
(few AU; Font et al. 2004; Liffman 2003) is being evaporated away faster than it can be 
replenished by viscous accretion. As a result, the inner disk is decoupled from the outer disk 
and accretes onto the star, leaving behind a true "inner hole" in the gas distribution and no 
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emission from the region within i?hoie- 

How do the properties of the gaseous inner disk of TW Hya compare with these predic- 
tions? UV fluorescent H 2 and CO fundamental emission, which are both believed to probe 
the inner ~ 1 AU of T Tauri disks (see Najita et al. 2007b for a review), have been detected 
from the inner disk of TW Hya (Herczeg et al. 2002; Rettig et al. 2004; Salyk et al. 2007; 
Najita et al. 2007b). The spectroastrometric study of Pontoppidan et al. (2008) further 
demonstrates that the CO emission arises from a rotating disk close to the star (~ 0.1 AU). 
Thus, the inner disk is not completely cleared of gas, a result that is consistent with the 
ongoing stellar accretion in the system and inconsistent with the EUV-driven photoevapora- 
tion scenario. (The situation is somewhat different if photoevaporation is driven by X-rays. 
Recent work by Ercolano and collaborators finds that photoevaporation driven by X-rays 
rather than EUV can create a transition-like SED at much higher disk masses and accretion 
rates than in the EUV-driven photoevaporation case. Under these conditions, the star may 
continue to accrete at a measurable rate for a longer fraction of the system lifetime after 
disk clearing has begun; Owen et al. 2010.) 

The Spitzer spectrum allows us to probe the gaseous disk at radii beyond the region 
probed by UV fluorescent H 2 and CO fundamental emission, because the features in the SH 
spectrum probe emission from cooler gas. In the spectrum of AA Tau, a typical T Tauri 
star that possesses an optically thick inner disk, the molecular features detected in SH (H2O, 
OH, C2H 2 , HCN, CO2), probe gas with temperatures of 300-600 K and disk emitting areas 
corresponding to a few AU in radius, i.e., the inner planet formation region of the disk (Carr 
k Najita 2008). 

One of the most striking aspects of the Spitzer spectrum of TW Hya is the weak molecu- 
lar emission compared to that seen in classical T Tauri stars such as AA Tau (Carr & Najita 
2008; Salyk et al. 2008; Carr & Najita, in preparation; Pontoppidan et al., in preparation). 
The scenario of a gap carved in the gaseous disk by an orbiting giant planet would appear to 
naturally lead to the suppression of the molecular emission from the inner disk region, like 
that observed. With the size of the optically thin region in the TW Hya system (~ 4AU) 
suggesting that an orbiting planet would reside at a few AU, we might expect the inner few 
AU region of the disk to be mostly cleared of gas. 

In comparison, because it suggests a continuous gaseous disk with no continuum opacity 
in the inner region, the grain growth scenario might predict similar or stronger mid-infrared 
molecular emission than would be observed in a continuous disk of gas and dust such as 
AA Tau, as long as the molecular abundances of the inner disk are not significantly altered 
by the optically thin nature of the inner disk. Photodissociation may alter the molecular 
abundances in the disk, although molecular shielding (e.g., by CO, H 2 , OH, and H 2 0) may 
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limit its impact. Another caveat is that molecules may be abundant in the inner disk of 
TW Hya, but molecular emission might not be excited to the level seen in AA Tau because 
of the lower accretion rate of TW Hya compared to typical T Tauri stars. A low accretion 
rate may imply a reduced rate of mechanical heating in the disk atmosphere. A low rate 
of mechanical heating can reduce the column density of warm water in the disk atmosphere 
(Glassgold et al. 2009) and lead to weak water emission. 

There is observational motivation to consider these issues. Photodissocation by stellar 
UV photons may play a role in reducing molecular abundances in the disk, and the hot OH 
we see may provide evidence for the photodissociation of water in this system (sections 3.5 
and 4.2). The lack of grains in a gas-rich inner disk may also limit the molecular abundances 
in this region. Glassgold et al. (2009) have shown that processes such as the formation of H 2 
on warm grain surfaces and mechanical heating associated with accretion can significantly 
enhance the abundance of water in the warm atmospheres of inner disks up to the levels 
observed in classical T Tauri disks. While it is yet unclear whether either of these is the 
dominant or critical factor in explaining the observed water abundances in classical T Tauri 
disk atmospheres, if H 2 formation on grains is an important requirement for the formation 
of H 2 in the warm disk atmosphere, the absence of grains in the inner disk region of a 
transition object may limit the H 2 emission from this region. 

Thus, the lack of strong MIR molecular emission from TW Hya indicates that the inner 
planet formation region of its gaseous disk (< 5AU) has evolved away from the classical T 
Tauri norm. But further work is needed to determine whether that evolution is the result of 
the formation of a giant planet or other physical or chemical processes. 

4.2. Origin of the Detected Emission 

4-2.1. Molecular Emission 

While strong molecular emission is not detected from TW Hya, we do detect a rich 
spectrum of weaker emission lines. Unlike the Spitzer spectrum of AA Tau (Carr & Najita 
2008), which shows a rich spectrum of H 2 0, OH, HCN, C 2 H 2 , and C0 2 emission, and the 
strong H 2 emission that is detected from the classical T Tauri stars DR Tau and AS 205A 
(Salyk et al. 2008), the molecular emission we detect from TW Hya in the 10-20 /zm region 
is dominated by radicals (OH, and possibly CH 3 ) and ions (HCO + ). 

Is the difference in the species detected in TW Hya vs. classical T Tauri stars due to 
the physical truncation of the disk, as discussed in the previous section? Alternatively, or 
in addition, does the difference result from the higher X-ray and UV irradiation field of TW 
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Hya? TW Hya has a nearly unique X-ray spectrum (Drake 2005) that shows unusually bright 
soft X-ray emission compared to other T Tauri stars (Robrade & Schmitt 2006). TW Hya 
is also known as a bright FUV emission source (Herczeg et al. 2004). Bergin et al. (2004) 
quote an FUV flux (including Ly a) of Go = 3400 at 100 AU for TW Hya, in comparison 
with values of 240, 340, and 1500 for DM Tau, GM Aur, and LkCal5, respectively. One 
hypothesis is that the stronger UV and soft X-ray field of TW Hya results from its lower 
accretion rate (lO~ 9 M yr _1 ; Muzerolle et al. 2000; Alencar & Basri 2000; Herczeg et al. 
2004), which unburies the accretion shock and allows more energetic photons to emerge 
from the "sides" of the accretion column (Drake 2005; Ardila 2007; Ardila & Johns-Krull 
2009). 

Perhaps this higher intensity photon field is responsible for some aspects of the emission 
that we see. A stronger UV field from TW Hya may dissociate H 2 in favor of OH. The hot 
OH emission that we detect may provide evidence for this process in action. This situation 
may therefore be similar to that of the Herbig Ae stars studied by Mandell et al. (2008), who 
suggested that the presence of OH emission and lack of H2O emission from these sources 
could be the result of photodissociation by stellar UV. 

If it arises from the photodissociation of water, the hot OH emission we observe may help 
to explain the origin of the high OH column density that is inferred for the disk atmospheres 
of classical T Tauri stars. Carr & Najita (2008) inferred a (line-of-sight) OH column density 
of 8 x 10 16 cm -2 at a temperature of ~ 500 K based on their analysis of the Spitzer spectrum 
of AA Tau. For a system inclination of i = 75 for AA Tau (Bouvier et al. 2007), the measured 
OH column density corresponds to a vertical column density of 2 x 10 16 cm~ 2 . Even larger 
OH column densities of 2 x 10 17 cm~ 2 were inferred for two classical T Tauri stars based on 
3/zm spectra (Salyk et al. 2008). 

In comparison, Glassgold et al. (2009) find that X-ray irradiated disks will produce warm 
(300-1000 K) OH (vertical) column densities of 1 x 10 14 — 8 x 10 14 cm~ 2 , with the specific 
value depending on the details of the model. Other chemical calculations for disks also fail to 
produce the large OH column densities (Agundez et al. 2008) or high OH/H 2 abundances 
(Woods & Willacy 2009) that are observed. Photodissociation of water to produce OH is 
one possible explanation for the much larger column density of OH that is observed (Bethell 
& Bergin 2009; Glassgold et al. 2009). The hot OH emission that we observe may provide 
observational support for that perspective. 

When water in the disk absorbs dissociating UV photons, the UV photons can both 
produce hot (prompt) OH emission as well as heat the disk atmosphere. After it has emitted 
the prompt emission, the OH produced in this process may relax to a thermal distribution 
(at the gas temperature) and continue to emit thermally. A UV irradiated disk may thereby 
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produce both thermal and prompt emission, as seen in a source such as AA Tau (section 
3.5). Detailed modeling is needed to determine if the OH emission observed in a system like 
AA Tau is consistent in detail with this picture. 

Where would such water be located in the TW Hya system, given that no water emission 
is observed at 10-20 pm? One possibility is that it resides in the optically thin region of 
the disk, but it is not excited enough or warm enough to emit in the 10-20 /im region. 
Alternatively, the water may be present in the outer, optically thick disk where it is too cool 
to produce much emission at 10-20 /im. To address this issue, it would be interesting to 
measure velocity resolved OH line profiles in order to determine if the OH emission arises in 
a disk and at what radii. 

The irradiation environment may also play a role in the HCO + emission that we detect. 
In their earlier study of organic molecules in the outer disk of TW Hya, Thi et al. (2004) 
found the abundance ratio of HCO + /CO to be ~ 10 times larger in TW Hya than in two 
other transition disks (LkCa 15 and DM Tau). Thi et al. suggested that the high ratio of 
HCO + /CO may result, in part, from the high X-ray flux of TW Hya. The high X-ray flux 
of TW Hya may also contribute to the prominence of HCO + emission in the mid-infrared. 

4-2.2. Atomic Emission 

As discussed in section 3.4, the [Nell] and [Nelll] fluxes we measure agree well with 
predictions for gaseous disks irradiated by stellar X-rays (see section 3.4). EUV irradiation 
of the disk may also contribute to some extent, although a soft UV spectrum is needed to 
be consistent with the large [NeII]/[NeIII] ratio observed. This result is consistent with the 
picture of a disk origin for the (spatially compact) [Nell] emission from T Tauri stars that is 
probed by high resolution spectroscopy. As discussed in Najita et al. (2009) based on high 
resolution spectroscopy of [Nell] emission from AA Tau and GM Aur, and Herczeg et al. 
(2007) in the context of TW Hya specifically, the [Nell] emission line is symmetric, centered 
at the radial velocity of the star, and has line widths that are plausibly explained by disk 
rotation. The possibility of an origin in jets or outflows, illustrated in the case of the T Tau 
triplet (van Boekel et al. 2009), is less likely here, because TW Hya is not a known outflow 
source. 

This perspective is basically confirmed, although with a twist, in the more recent study 
of [Nell] in TW Hya using VISIR at the VLT. Pascucci & Sterzik (2009) interpret their 
measurement of the [Nell] line profile of TW Hya, which shows a Gkms" 1 blueshift in the 
emission centroid, as indicating an origin in a low-velocity photoevaporative flow from the 
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disk. This result is related to the idea of a disk origin for [Nell] emission in that the emission 
still arises from disk gas, although it is gas that is gravitationally unbound rather than in 
hydrostatic equilibrium. 

The HI emission we detect is not as well known or understood. While the 12.37 jum 
HI 7-6 line has been reported previously from T Tauri stars (e.g., Pascucci et al. 2007; see 
also Ratzka et al. 2007 for TW Hya), the weaker emission lines have not been reported 
previously. Perhaps such emission lines are commonly present in T Tauri stars, but they 
are obscured by stronger molecular emission and the weakness of the molecular emission in 
TW Hya makes the weaker HI lines easier to detect in this source. On the other hand, TW 
Hya is an energetic UV and X-ray emission source, as discussed above, and the emission we 
detect may be more unusual as a result. 

What component of the T Tauri system could potentially produce the HI emission 
that we observe? Ratzka et al. (2007) suggested that the HI lines most likely originate 
from either an accretion shock or the stellar corona, although they were not specific on the 
details. Pascucci et al. (2007) have argued that the HI 6-5 and HI 7-6 emission detected from 
other low accretion rate T Tauri stars are are too strong to be explained by magnetospheric 
accretion flows. 

Following the lead of previous studies that interpreted the Br7 and Balmer lines of 
T Tauri stars as arising in a magnetospheric accretion flow (Calvet & Hartmann 1992; 
Hartmann, Hewett, & Calvet 1994; Najita et al. 1996; Muzerolle et al. 1998a,b; Folha & 
Emerson 2001), Bary et al. (2008) interpreted the line ratios of the high-lying Paschen and 
Brackett lines they observed at low spectral resolution as arising in a magnetosphere with 
unusually low temperature. As described by Bary et al., the electron densities predicted 
theoretically for funnel flows (Martin 1996; Muzerolle et al. 2001) are consistent with their 
line ratios. They did not discuss whether magnetospheres can reproduce the fluxes of the 
Paschen and Brackett lines they observe. 

The MIR HI lines impose further constraints on the origin of the HI emission. It is 
unclear whether the HI emission we observe can be produced by the accretion shocks that 
are discussed in the literature. The accretion shock in the TW Hya system is expected to have 
a height of ~ 1-1000 km for the post-shock region (Ardila & Johns-Krull 2009; Ardila 2007) 
and a density of ~ 10 13 cm" 3 (Drake 2005; Ardila 2007). The pre-shock region in typical 
T Tauri stars is characterized as having a temperature of ~ 10,000 K, a height of ~ 1000 
km, an average density of 10 11 — 10 15 cm" 3 , and a filling factor of / ~ 10~ 3 — 10~ 2 (Calvet 
& Gullbring 1998). These densities are high enough to produce inverted HI populations 
(Hummer & Storey 1987). However, significant laser amplification may not occur or may be 
difficult to detect because of laser saturation and/or beaming effects (Strelnitzki et al. 1996; 
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Smith et al. 1997). Further work is needed to determine whether such high density, compact 
regions can produce bright HI emission with roughly case B line ratios. 

The line profiles of transitions connecting high lying HI levels may also provide some 
insights. One example is the Pa7 (6-3) profile of TW Hya (Edwards et al. 2006), which is 
superficially similar to a magnetospheric accretion profile: the profile is steeper on the red 
side than the blue and the blue wing extends to ~ 300 km s -1 . However, as noted by Edwards 
et al., the profile is not consistent in detail with predictions for magnetospheric accretion, 
e.g., the profile is narrower and more centered than predicted profiles (e.g., Muzerolle et 
al. 2001). Such narrow Pa7 profiles are found commonly for low accretion rate systems 
(Edwards et al. 2006). 

Edwards et al. note a similarity with the narrow component of metallic emission lines 
(e.g., Fel and Fell), which are also the dominant component of the profile in low accretion 
sources such as TW Hya. The narrow and broad components are interpreted as arising 
in different physical regions, with the broad component having a significant contribution 
from funnel flows and an additional physically extended component from winds. The narrow 
component is associated with an accretion shock at the footpoint of the funnel flow (Edwards 
et al. 2006), which may not be the region from which the MIR HI emission originates, as 
noted above. 

Disk atmospheres and/or photoevaporative flows may also contribute detectable HI 
emission. The disk atmosphere and photoevaporative flow that results from the X-ray ir- 
radiation of the disk, as described by Ercolano et al. (2008), produces a line luminosity of 
4 x 10 25 ergs~ 1 in the HI 7-6 (12.37 /zm) line. In an updated calculation, Ercolano et al. 
(2009) solve for the hydrostatic structure of the disk and include irradiation by both EUV 
and X-rays. In this case, the HI 7-6 flux is revised upwards to 4 x 10 26 ergs _1 in FSOHOLxl 
(B. Ercolano, private communication 2009), a model that is fairly optimistic in assuming 
no attenuation of the EUV flux by foreground HI absorption (e.g., as might be present in a 
magnetosphere) . We observe a yet larger luminosity of 1.8 x 10 28 ergs _1 , the average of the 
flux in SHI and SH2. 

A more general question is whether the neon and HI emission can arise from the same 
physical region of the TW Hya system. The picture of [Nell] and [Nelll] emission arising in 
a mostly neutral region (e.g., in an X-ray irradiated disk atmosphere; section 3.4) contrasts, 
at least superficially, with the nominal picture of HI recombination line emission. Case B 
hydrogen line ratios are typically interpreted as indicating formation in a fully ionized region. 
(While a high ionization fraction is clearly not required to produce recombination emission, 
it helps to maximize the emission from a given mass of gas.) 
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The mismatch in these conditions seems to show up in other models of EUV and X-ray 
irradiated disks as well. Hollenbach & Gorti (2009) show that although irradiation by soft 
EUV can account for the flux of the [Nell] emission observed from T Tauri stars, it cannot 
produce the observed HI 7-6 line luminosity: the predicted HI 7-6 flux is < 1% of the 
[Nell] flux. As a result, Hollenbach & Gorti (2009) argue that the HI emission arises from a 
different region than the [Nell]. As they note, a high density region is one possibility since 
it may produce significant HI emission without producing additional (unwanted) [Nell]. 
Considerations of HI population inversions and laser amplification, and whether the HI 
emission would have (optically thin) case B line ratios, as described above, may limit the 
allowable range of densities. 

Given the current uncertainty regarding the possible origin of the HI emission, it ap- 
pears that high resolution spectroscopy of the MIR HI lines would likely provide valuable, 
and much needed, insights into the origin of the emission. Since the 11.31 /iin and 12.37/zm 
lines are comparable in brightness to the [Nell] line, such a study should be feasible cur- 
rently. If the mid-infrared lines are found to be as broad as the Pa7 line, the role of high 
velocity phenomena such as magnetospheric infall would be indicated. If the lines are nar- 
row (> 30kms _1 ) and/or with a small blueshift (~ lOkms^ 1 ), disk atmospheres and/or 
photoevaporative flows would be indicated. 

Both components might in fact contribute. As noted above, the Pa7 profiles of Ed- 
wards et al. (2006) are narrower and more centered in velocity than magnetospheric emis- 
sion models. A combination of a high- velocity magnetospheric component and a low velocity 
disk/photoevaporative component may better account for the observed profiles than a mag- 
netosphere alone. 

5. Summary and Future Directions 

The Spitzer spectrum of TW Hya is strikingly different from that of other classical T 
Tauri stars reported in the literature, showing weak molecular emission in the 10 — 20 /iin 
region and displaying no strong emission features of H2O, C2H2, or HCN. The difference 
indicates that the inner planet formation region (< 5AU) of its gaseous disk has evolved 
away from the classical T Tauri norm. With the SED of TW Hya indicating that the dust 
disk within ~ 4 AU of the star has become optically thin in the continuum, possibly because 
it has been cleared by an orbiting giant planet (Calvet et al. 2002), it is tempting to infer 
that TW Hya shows weak molecular emission in the mid-infrared because the planet has 
created a gap in the gaseous disk as well. We discussed in section 4 some of the issues that 
need to be investigated in order to determine whether the difference we observe is the result 
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of giant planet formation or other processes. 

Given our limited understanding of the factors that govern the emission spectra of T 
Tauri disks, it may be useful to take an empirical approach in exploring whether the lack 
of molecular emission seen from TW Hya is a consequence of the physical evolution of the 
disk, its chemical evolution, or an excitation effect. For example, to explore the possibility 
of an excitation effect, we can compare the Spitzer spectrum of TW Hya with those of other 
non-transition objects with comparable accretion rates. If non-transition objects with low 
accretion rates also lack strong molecular emission, the spectrum of TW Hya would not be 
unusual for its accretion rate and the lack of molecular emission may be the result of poor 
excitation. However, if strong molecular emission is observed in other low accretion rate 
systems, the absence of such emission in TW Hya would suggest a gap in its gaseous disk or 
a possible chemical effect. We will take this approach in a future study. 

Similarly, it would be interesting to explore whether other transition objects also show 
weak molecular emission compared to classical T Tauri stars. If transition objects with 
much weaker UV fluxes than that of TW Hya (e.g., DM Tau, GM Aur) also show a deficit 
of molecular emission, that would suggest a more dominant role for SED evolution (a deficit 
of grains; possible clearing by a giant planet), rather than photochemistry, in accounting for 
the difference in the spectra. We will report on this in future publications. 

While strong molecular emission is not detected from TW Hya, we do detect a rich 
spectrum of emission lines of atoms (HI, [Nell], and [Nelll]) and molecules (H 2 , OH, C0 2 , 
HCO + , and possibly CH 3 ). One of the most intriguing is the OH emission, which is hot and 
may result from the UV photodissociation of water. A more detailed analysis of the OH 
emission spectrum may be able to determine whether it is produced by photodissociation. 
The properties of the molecular emission from TW Hya, both the OH and other molecules, 
will be analyzed in greater detail in a future study. 

Because we detect multiple HI lines, we can show that the HI emission from TW Hya 
has a recombination spectrum. In contrast to the neon emission from TW Hya, which can 
be well accounted for by (primarily) stellar X-ray irradiation of the disk, the physical origin 
of the HI emission is difficult to identify. As discussed in section 4, magnetospheres, disk 
atmospheres and/or photoevaporative flows could plausibly contribute to the emission and 
multiple components may play a role. High resolution spectroscopy of the brightest HI lines 
(HI 7-6 and HI 9-7) would likely provide valuable insights into the origin of the emission. 

We are grateful to Al Glassgold, Barbara Ercolano, and David Ardila for interesting 
and useful discussions regarding the interpretation of the observations. JN thanks for their 
generous hospitality Tom Soifer and the Spitzer Science Center, where much of the analysis 
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Table 1. Observing Log for TW Hya SH Datasets 



Dataset 


Date 


AOR 


Ramp 


On-Source 


Off-Source 








(sec) 


(sec) 


(sec) 


GTO 


2004-01-04 


3571456 


6 


12 





SHI 


2006-07-01 


18017792* 


30 


120 


120 


SH2 


2008-01-07 


24402944 


30 


600 


300 



*AOR for sky observations: 18018048, 18018304. 
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Fig. 1. — Spitzer IRS short-high spectra of TW Hya taken at three epochs: SHI (top), SH2 
(bottom), and the GTO spectrum (middle). The GTO spectrum is offset by +0.15 Jy and 
SHI by +0.25 Jy for clarity. 
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Fig. 2. — Average of the SHI and SH2 spectra of TW Hya from which a polynomial fit to the 
continuum has been subtracted. Emission features of atomic lines (HI, [Nell], and possibly 
[Nelll]) and molecules (H 2 , OH, C0 2 , HC0+, and possibly CH 3 ) are detected. The [NelH] 
and CH 3 features are clearer in the SH2 and SHI spectra, respectively (see Figure 3). 
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Fig. 3. — The region around the [Nelll] line in the SH2 spectrum (top) and the CH 3 emission 
band in the SHI spectrum (bottom). A representative errorbar (± 1-a) is shown in the upper 
left corner of each panel. 




Fig. 4. — FWHM of the Gaussian fits to OH features (large red dots) and [Nell] and un- 
blended HI lines (black squares) in the SH2 spectrum. Where errorbars are not apparent, 
they are smaller than the size of the symbol. The nominal resolution of IRS in the short- 
high mode (R = 750; black line) is shown for comparison. While the HI and [Nell] lines 
are spectrally unresolved, the widths of the OH features grow with wavelength due to the 
increasing spread in the wavelengths of individual OH lines that make up a given feature. 
This is illustrated by the agreement with the small blue dots, which show this spread in 
wavelengths added in quadrature with the spectral resolution at the central wavelength of 
each feature. 
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Fig. 5. — Ratio of line fluxes (top), equivalent widths (middle), and Gaussian FWHM 
(bottom) in SHI compared to SH2 for HI (black dots), [Nell] (red dot), H 2 (open blue 
squares), OH (solid blue triangles), and HCO + and CH 3 (open blue triangles). The fluxes 
of the HI lines are larger by ~ 2 in flux in SHI compared to SH2. The equivalent width of 
the [Nell] line is smaller in SHI than SH2. 
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Fig. 6. — HI line luminosities from the SHI spectrum (dots with error bars) compared with 
relative fluxes from case B recombination (triangles). The HI 15-9, HI 15-10, and HI 14-10 
lines are not detected, consistent with their low anticipated fluxes. 
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Fig. 7. — HI line luminosities from the SH2 spectrum (dots with error bars) compared with 
relative fluxes from case B recombination (triangles). The HI 8-7 line is weaker at this epoch 
and not in as good agreement with the case B relative fluxes. The HI 15-9, HI 14-9, HI 
15-10, and HI 14-10 lines are not detected, consistent with their low anticipated fluxes. 
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Table 2. 


Continuum and Errors 










CTO 




SHI 




SH9 




Feature 


Wavelength 


Continuum 


Error 


Continuum 


Error 


Continuum 


Error 




(/im) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


HI 9-7 


11.31 


0.68 


0.013 


0.81 


0.0037 


0.62 


0.0019 


Nell 


12.81 


0.51 


0.011 


0.61 


0.0030 


0.47 


0.0017 


co 2 


14.97 


0.76 


0.015 


0.89 


0.0036 


0.69 


0.0019 


H2 S(l) 


17.03 


1.28 


0.023 


1.47 


0.0063 


1.18 


0.0050 


HI 8-7 


19.06 


1.85 


0.029 


2.09 


0.012 


1.70 


0.0084 



Table 3. OH Features Detected 



A/ /im v Transitions 



10 


.81 





3/2 


R 


29. 


,5, 


1/2 


R 


28. 


,5 


11 


,06 





3/2 


R 


28. 


,5, 


1/2 


R 


27. 


.5 


11 


.31 





3/2 


R 


27. 


.5, 


1/2 


R 


26. 


.5 


11 


.61 





3/2 


R 


26. 


.5. 


1/2 


R 


25. 


,5 


11 


.92 





3/2 


R 


25. 


.5, 


1/2 


R 


24. 


,5 


12. 


.65 





3/2 


R 


23. 


.5, 


1/2 


R 


22. 


.5 


13 


.08 





3/2 


R 


22. 


,5, 


1/2 


R 


21. 


.5 


13 


.54 





3/2 


R 


21. 


.5, 


1/2 


R 


20. 


,5 


14. 


.08 





3/2 


R 


20. 


.5, 


1/2 


R 


19. 


,5 


14. 


,65 





3/2 


R 


19. 


.5, 


1/2 


R 


18. 


.5 






1 


3/2 


R 


20. 


.5, 


1/2 


R 


19. 


,5 


15. 


.30 





3/2 


R 


18. 


,5, 


1/2 


R 


17. 


,5 






1 


3/2 


R 


19. 


,5, 


1/2 


R 


18. 


,5 


16 


.02 





3/2 


R 


17. 


,5, 


1/2 


R 


16. 


.5 


16 


.84* 





3/2 


R 


16. 


.5, 


1/2 


R 


15. 


.5 


17. 


.78 





3/2 


R 


15. 


,5, 


1/2 


R 


14. 


.5 



*Blended with the HI 12-9 line. 
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Table 4. Emission Line Properties in the SHI Spectrum 



feature 


A 


r lux 


err 


r WHM 


+err 


—err 


EW 


err 




(H 


( 1 n — 14 

(10 erg 


-1 — 2\ 

;s cm J 


(10 /im) 






(10" 3 /im) 


HI 12-8 


10.502 


1.07 


0.54 


2.2 


2.1 


0.9 


0.47 


0.24 


OH 


10.809 


0.82 


0.51 


2.8 


3.5 


1.9 


0.38 


0.24 


OH 


11.077 


1.26 


0.73 


5.1 


5.2 


2.4 


0.62 


0.36 


HI 9-7 


11.310 


3.74 


0.38 


1.6 


0.2 


0.1 


1.96 


0.20 


OH 


11.609 


0.69 


0.40 


2.0 


1.9 


1.2 


0.42 


0.24 


OH 


11.922 


0.66 


0.30 


1.8 


1.6 


0.7 


0.46 


0.21 


HCO + 


12.064 


1.51 


0.41 


2.8 


1.1 


0.6 


1.13 


0.31 


HoSf2)+OH 

-L -L^ kJ I 1 | Vy J. J. 


12.276 


1.49 


0.42 


3.4 


1.4 


0.9 


1.22 


0.34 


HI 7-6 


12.374 


8.10 


0.33 


2.1 


0.9 


1.2 


6.75 


0.28 


HI 14-9 


12.587 


0.51 


0.25 


1.5 


1.6 


0.7 


0.43 


0.21 


OH 


12.663 


0.37 


0.28 


1.3 


2.0 


0.7 


0.33 


0.25 


Nell 


12.814 


4.95 


0.29 


1.8 


0.1 


0.1 


4.43 


0.26 


OH 


13.069 


0.26 


0.19 


1.3 


2.4 


0.9 


0.23 


0.18 


OH 


13.541 


0.80 


0.33 


3.4 


1.9 


0.9 


0.73 


0.30 


OH 


14.074 


0.31 


0.23 


2.3 


4.7 


1.9 


0.28 


0.20 


HI 13-9 


14.183 


0.89 


0.30 


3.1 


1.6 


1.2 


0.79 


0.27 


OH 


14.645 


1.23 


0.42 


4.7 


2.1 


1.4 


1.06 


0.36 


co 2 


14.974 


0.70 


0.47 


5.4 


6.1 


3.5 


0.59 


0.39 


OH 


15.302 


0.43 


0.26 


1.7 


2.4 


0.9 


0.34 


0.21 


OH 


16.013 


1.56 


0.59 


6.4 


4.2 


1.6 


1.12 


0.42 


HI 10-8 


16.215 


1.78 


0.51 


3.8 


1.4 


1.2 


1.23 


0.36 


CH 3 


16.482 


0.45 


0.33 


1.9 


4.7 


1.6 


0.30 


0.23 


HI12-9+OH 


16.859 


2.53 


0.74 


8.0 


3.5 


2.1 


1.69 


0.49 


H 2 S(l) 


17.036 


1.21 


0.41 


2.0 


1.2 


0.5 


0.79 


0.27 


OH 


17.762 


1.38 


0.57 


4.7 


2.6 


2.1 


0.84 


0.35 


HI 8-7 


19.064 


3.06 


0.83 


3.1 


0.9 


0.9 


1.77 


0.48 
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Table 5. Emission Line Properties in the SH2 Spectrum 



Feature 


A 


Flux 


err 


FWHM 


+err 


—err 


EW 


err 




(H 


(10~ 14 erg: 


3" 1 cm" 2 ) 


(10- 2 /im; 


) 




(io- 3 


fim) 


HI 12-8 


10.508 


0.45 


0.24 


1.8 


1.9 


0.9 


0.26 


0.14 


OH 


10.8146 


0.38 


0.26 


2.7 


4.7 


1.4 


0.23 


0.16 


OH 


11.063 


0.43 


0.20 


1.5 


1.2 


0.7 


0.27 


0.12 


HI 9-7 


11.310 


1.60 


0.19 


1.4 


0.2 


0.2 


1.10 


0.13 


OH 


11.614 


0.15 


0.13 


0.9 


1.6 


0.8 


0.12 


0.11 


OH 


11.922 


0.52 


0.23 


2.9 


2.4 


1.2 


0.47 


0.21 


HCO+ 


12.064 


1.32 


0.21 


2.6 


0.5 


0.6 


1.30 


0.20 


H 2 S(2)+OH 


12.276 


1.14 


0.20 


2.9 


0.7 


0.7 


1.21 


0.21 


HI 7-6 


12.373 


4.53 


0.20 


2.1 


0.1 


0.1 


4.92 


0.21 


OH 


12.659 


0.45 


0.19 


2.2 


1.4 


0.7 


0.52 


0.21 


Nell 


12.813 


5.56 


0.16 


1.72 


0.0 


0.1 


6.47 


0.19 


OH 


13.081 


0.55 


0.21 


3.4 


1.9 


1.2 


0.65 


0.25 


OH 


13.547 


0.48 


0.20 


2.9 


1.9 


1.2 


0.56 


0.23 


Url 


1 A D7D 
14. U 1 U 


U.OZ 


U.zo 


4.z 


z.l 




U. /z 


U.ZO 


HI 13-9 


14.181 


0.20 


0.18 


2.6 


4.7 


1.6 


0.23 


0.21 


OH 


14.651 


0.94 


0.22 


4.5 


1.2 


1.4 


1.04 


0.25 


co 2 


14.972 


0.62 


0.17 


2.7 


1.1 


0.6 


0.67 


0.19 


OH 


15.294 


0.75 


0.28 


5.4 


2.8 


1.9 


0.77 


0.29 


Nelll 


15.558 


0.25 


0.19 


1.6 


2.1 


1.2 


0.24 


0.18 


OH 


16.016 


1.01 


0.35 


5.2 


2.6 


1.4 


0.92 


0.32 


HI 10-8 


16.212 


0.88 


0.26 


2.6 


1.2 


0.7 


0.77 


0.23 


HI12-9+OH 


16.861 


1.69 


0.58 


8.0 


3.8 


2.1 


1.41 


0.48 


H 2 S(l) 


17.035 


1.26 


0.31 


1.9 


0.7 


0.5 


1.03 


0.25 


OH 


17.767 


1.08 


0.57 


6.2 


5.9 


2.6 


0.82 


0.43 


HI 8-7 


19.061 


0.71 


0.50 


2.1 


2.1 


1.4 


0.50 


0.35 
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Table 6. Emission Line Properties in the GTO Spectrum 



Feature 


A 


Flux 


err 


EW err 




(H 


(10~ 14 ergs 


_1 cm~ 2 ) 


(10- 3 /im) 


HI 9-7 


11.306 


2.92 


1.53 


1.83 0.96 


HCO+ 


12.060 


2.39 


1.44 


2.12 1.28 


H 2 S(2) 


12.272 


1.55 


1.00 


1.51 0.97 


HI 7-6 


12.372 


7.13 


1.37 


7.13 1.37 


Nell 


12.812 


5.90 


0.91 


6.29 0.97 



